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Density-functional theory is used to study the geometric and electronic structure 
of cationic Si]^g clusters with a Ti, V or Cr dopant atom. Through unbiased global 
geometry optimization based on the basin-hopping approach we confirm that a Frank- 
Kasper polyhedron with the metal atom at the center represents the ground-state 
isomer for all three systems. The endohedral cage geometry is thus stabilized even 
though only VSi^g achieves electronic shell closure within the prevalent spherical 
potential model. Our analysis of the electronic structure traces this diminished role 
of shell closure for the stabilization back to the adaptive capability of the metal- 
Si bonding, which is more the result of a complex hybridization than the orginally 
proposed mere formal charge transfer. The resulting flexibility of the metal-Si bond 
can help to stabilize also "non- magic" cage-dopant combinations, which suggests 
that a wider range of materials may eventually be cast into this useful geometry for 
cluster-assembled materials. 
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I. INTRODUCTION 



Doping with endohedral metal atoms appears as a remarkable avenue to tailor the in- 
trinsic properties of silicon clusters."^ In contrast to the compact geometries caused by the 
preferable sp^ bonding in pure Si clusters-"^, the incorporation of even a single impurity 
atom can lead to the stabilization of otherwise unfavorable cage-like structures'""^. As in 
clathrates^^ or carbon nanostructures, these cages then represent appealing symmetric and 
unreactive building blocks for novel cluster-assembled materials with engineered properties. 
A prerequisite to a systematic synthesis of such materials are simple rules that rationalize 
which metal dopants stabilize cages and of which size. While this has been controversially 
discussed, one commonly agreed criterion for highly stable so-called "magic" clusters is geo- 
metric and electronic shell closure .-liiii^ Here, the electronic manifold of a highly symmetric 
cage is viewed as states in a spherical potential, and particular stability is expected, if the 
electrons fill any one of the angular momentum shells, i.e. Is (2e~), Ip (6e~), Id (10e~), 1/ 
(14e~), 2s (2e-), Ig (ISe"), 2p (6e-), 2d (lOe") etc' 

For a 16 Si atom endohedral Frank-Kasper (FK) polyhedron'*' "magicity" would hence 
be predicted for a dopant atom donating 4 valence electrons, as the resulting 16 x 4 + 4 = 68 
electrons achieve closure of the 2d shell. Within this model the known high stability of VSi^'g 
is thus naturally explained, if the nature of the bonding in the cluster is viewed in terms 
of a full formal charge transfer, i.e. "VSi^g = V^^ + Si^^^".'" Recently, however, Lau 
et al. deduced from their X-ray absorption spectroscopy data that also TiSi^g and CrSi^^g 
with one valence electron less and more, respectively, stabilize in a cage geometry, with 
furthermore a highly similar local electronic structure around the dopant atom compared 
to the classic VSi^'g system.— These findings motivate the present theoretical study, which 
follows a twofold goal. First, we perform a first-principles global geometry optimization 
of the three cluster systems MSi^^g (M =Ti,V,Cr) to firmly establish that the endohedral 
FK cage indeed represents the ground state geometry for all three dopant atoms. Second, 
we analyze the obtained electronic structure to obtain a more qualified view on the nature 
of the chemical bonding and elucidate the mechanism that stabilizes the cage despite the 
differing number of valence electrons in the three systems. 
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II. THEORY 



All calculations have been performed with the all-electron full-potential density-functional 
theory (DFT) code FHI-aims^^. Electronic exchange and correlation was treated within the 
generalized-gradient approximation functional due to Perdew, Burke and Ernzerhof (PBE)^^. 
For comparison single-point calculations at the optimized PEE geometries were also per- 
formed on the hybrid functional level using the B3LYP^^ and PBEG^'* functionals, without 
obtaining results that would lead to any conclusions different to the ones derived and pre- 
sented below on the basis of the PEE data. FHI-aims employs basis sets consisting of 
atom-centered numerical orbitals. All sampling calculations are done with the "tier2" basis 
set, which contains 43 basis functions for Si, 67 basis functions for Ti, 88 basis functions 
for V, and 88 basis functions for Cr, respectively. The numerical integrations have been 
performed with the "tight" settings, which correspond to integration grids with 85, 97, 99, 
and 101 radial shells for Si, Ti, V, and Cr, respectively, in which the number of integration 
points is successively decreased from 434 for the outermost shell to 50 for the innermost 
one.— For the ensuing electronic structure analysis of the optimized geometries the electron 
density was recomputed with an enlarged "tier3" basis set, which contains 64 basis functions 
for Si, 103 basis functions for Ti, 115 basis functions for V and 124 basis functions for Cr. 
Systematic convergence tests indicate that these settings are fully converged with respect to 
the target quantities (energetic difference of different isomers in the sampling runs; electron 
density distribution in the electronic structure analysis). This holds in particular for a cen- 
tral quantity of our analysis, the radial electron density distribution of the different doped 
cages. This quantity is defined as 



where n{r, 9, (p) is the electron density at a given point at spherical coordinates (r, 6, 0) away 
from the central dopant atom at r = 0. To build a radial distribution of the electron density 
we calculate the surface integral, eq. (1), for a set of spheres of increasing radii, and then 
plot the obtained values as a function of the sphere radius. The numerical integration is 
hereby performed using a cubic (400 x 400 x 400) volumetric data grid with 0.02 A voxel 
width. The chosen finite integration radius and angle steps equal 0.02 Bohr and 7r/360, 
respectively. 




(1) 
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Local structure optimization is done using the Broyden-Fletcher-Goldfarb-Shanno method^^ 
relaxing all force components to smaller than 10~^eV/A. To make sure that the cage-like 
geometry indeed represents the ground-state structure for all three dopant atoms we relied 
on basin-hopping (BH) based global geometry optimizatioii^^"— . Within the BH idea the 
configuration space is explored by performing consecutive jumps from one local minimum 
of the potential energy surface (PES) to another. To achieve this, positions of atoms in 
the cluster are randomly perturbed in a so-called trial move followed by a local geometry 
optimization which brings the system again into a local PES minimum. A Metropolis-type 
acceptance rule is used to either accept or reject the jump into the PES minimum reached by 
the trial move. As specific BH implementation we chose collective as well as single-particle 
trial moves, in which all atoms (collective move) or a randomly picked atom (single-particle 
move) is displaced in a random direction. Two different starting points were used for all opti- 
mization runs: 1) All atoms are randomly positioned inside a box of dimension (9 x 9 x 9) A^, 
or 2) the solution for the Thomson-problen>22 (how to put point charges on a sphere with 
minimal energy) is employed to position the Si atoms and then the doping metal atom is 
added at the center. Typical BH runs comprised of the order of 100 accepted trial moves, 
and unanimously identified the cage geometry as lowest energy structure regardless of the 
specific settings employed for the Metropolis rule or the single-particle/collective moves. 

III. RESULTS 

A. Cage-like ground state geometry 

In contrast to our preceding work on Si2o fuUerenes^S^ our extended unbiased configuration 
searches confirm that the endohedral FK cage indeed represents the ground-state isomer for 
MSi^g with all three dopant atoms, cf. Fig. [1] The "non-magicity" in case of Ti and Cr 
doping only expresses itself in form of a much reduced energetic gap to the second lowest 
energy isomer identified in the BH runs: For VSi^g this gap amounts to LOOeV, whereas for 
TiSi^g and CrSi^g it is only 0.01 eV and 0.08 eV, respectively. Within a 1.00 eV range above 
the identified ground state we correspondingly found about 15 inequivalent isomers for the 
latter two systems. In the Cr-doped case all of them are capped CrSi^^g cages, for Ti more 
compact TiSi^'g cages are found within 0.1 eV above the ground state. Above this mostly 
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FIG. 1. Ball-and-stick views of the identified ground-state FK cage geometries. 

capped TiSi^5 structures are identified. 

The incomplete shell-closure of TiSi^g and CrSij/g also shows up in the symmetry of the 
FK cage. Whereas the "ideal" VSi^g cluster exhibits perfect symmetry, the cages with Ti 
and Cr dopants only exhibit Ci symmetry. The distortions away from perfect Td symmetry 
are, however, only minor, as can best be seen from the M-Si bond distances within the cage. 
In a perfect FK polyhedron these distances fall into two closely spaced shells: One with four 
Si atoms that form a perfect tetrahedron, and slightly beyond that another one with 12 Si 
atoms that are all equidistant from the encapsulated metal atom. In the VSi^g cluster, these 
two shells are located at distances of 2.54 A and 2.81 A, respectively. In the less symmetric 
TiSi^'g and CrSi^'g geometries the distortions lift the degeneracies of the two shells and we 
instead find M-Si distances spread over a range of 2.64 A to 2.86 A for Ti and over a range 
of 2.50 A to 3.35 A for Cr, respectively. Overall this leads in case of Ti doping to a slightly 
increased average cage radius of 2.78 A, compared to the average M-Si distance of 2.74 A in 
both VSifg and CrSi^g. 

Overall, the geometric differences between the three cages are thus rather small. We 
furthermore verified that these differences have only insignificant effects with respect to the 
discussion of the electronic structure of the cage presented in the following. Our analysis 
is therefore for all three doped cages, as well as the empty Si cage based on the symmetric 
geometry obtained for VSi^'g, i.e. in all cases the geometry was kept as in VSi^g and only 
the electron density was each time self-consistently computed. This procedure facilitates the 
qualitative discussion of the nature of the bonding and of the concomitant charge redistribu- 
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FIG. 2. Total density of states (DOS) and DOS projected on the metal dopant for A) TiSi^, B) 
VSi^g, and C) CrSi^g. Panel D) directly compares the metal-projected DOS for the three cases to 
illustrate the varying degree of metal-Si hybridization. The zero-reference for the energy scale is 
the vacuum level, and the labels given to the different groups of states follow the notation of the 
spherical potential model (see text). 

tion as it allows to directly subtract the electron densities obtained for the different dopants 
and, because of the higher symmetry, makes the presentation of radial electron distributions 
averaged over the solid angle more meaningful. 



B. Spherical potential model 

The prevalent model to rationalize the stability of doped Si cage geometries is the spherical 
potential modet^*^"— which has been discussed in detail for the "magic" VSi^g cluster 
by Torres, Fernandez and BalbasA^. As a first step in our attempt to qualify the chemical 
bonding and stability in the "non-magic" FK clusters doped with Ti and Cr we first briefly 
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recapitulate the essentials of this discussion. The spherical potential model exploits the near 
sphericity of the ideal FK polyhedron, which suggests to classify the electronic states in shells 
of a determined radial and angular momentum quantum number. The computed density of 
states (DOS) of VSi]*g shown in Fig. |2J3 demonstrates that the Kohn-Sham states indeed 
group into the expected sequence (Is, Ip, Id, If, 2s, Ig, 2p, 2d, Ih, . . .), with the 68 valence 
electrons exactly achieving closure of the 2d shell. Also more subtle features like the splitting 
into the different tetrahedral (Td) sub-groups are perfectly obeyed, i.e. the different shells are 
sub-divided as s{ai),p{t2), d{e + t2), /(a2 + i^i + ^2), 5'(ai + e + + t2), h{e + ti + 2t2). Bonding 
to the transition metal dopant is predominantly expected via the vr-type orbitals with one 
radial node (2s, 2p, 2d), with hybridization following an approximate /-selection rule, i.e. the 
dopant 3d and 4s valence orbitals mix with Si cage d and s 7r-orbitals, respectively. The 
metal-projected DOS contained in Fig. [2j3 proves that also this feature of the spherical 
potential model is fully reproduced by the actual computation. 



However, these features are not a specificity of the "magic" VSiJg cluster, but instead 
inherent properties of the near-spherical FK polyhedral shape. As apparent from Fig. |2JA. 
and C essentially the same groupings of the Kohn-Sham states are equally obtained for the 
other two dopants, i.e. also here the electronic manifold is well described within the spherical 
potential model. Exactly as expected from the differing number of valence electrons, the 
only difference is that electronic shell closure is not achieved. In TiSi^g (with 67 valence 
electrons) the highest energy state of the 2d shell is unoccupied, and in CrSi^g (with 69 
valence electrons) the lowest energy state of the Ih shell is occupied. This lifts many of 
the degeneracies within the different electronic shells, but the overall structure in terms of 
angular momentum shells is still preserved. Furthermore, as confirmed by our first-principles 
sampling calculations the endohedral FK polyhedron still represents the lowest-energy isomer 
for the "non-magic" TiSi^g and CrSi^g. Electronic shell closure might thus be a criterion 
for enhanced stability, as e.g. reflected by the abundance of VSi^g in the experimental mass 
spectra. However, it is not a necessary condition to stabilize the endohedral cage geometry 
per se. 
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FIG. 3. Radial distribution, cf. eq. (1), of the electron density difference, An(r) = nygj+ (r) — 
ngj4- (r) — nv5+ (r) (solid line), where ?T-vsi+. ('^) electron density of the doped cage, ngj4- (r) the 

density of the empty Si cage, and nv5+ (r) the density of the V cation. If the formal charge transfer 
picture was correct, An(r) should be essentially zero throughout. Note the average cage radius, 
i.e. the position of the Si atoms, is at 2.74 A. Additionally shown are other charge combinations 
of the two fragments (Siie + V"*", dashed line; Si^g + V, dotted line). 

C. Charge transfer vs. hybridization 

Insight into the weakened role of electronic shell closure can come from a more qualified 
discussion of the nature of the chemical bond within the doped clusters. The simplified 
picture connected with the "magicity" of VSi^g assumes a formal charge transfer of all 
V valence electrons to the Si cage manifold. This "formal" view is readily checked by 
evaluating the difference of the actually computed electron density of VSi^g with respect to 
a mere superposition of the electron densities of an empty Si^g cage and a cation. If the 
formal charge transfer picture was correct, then this electron density difference should be zero 
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FIG. 4. Same as Fig. [3l but for TiSi^g (left panel) and for CrSi^i^g (right panel). 



throughout. Figure [3] shows this difference in form of the radial electron density distribution, 
i.e. averaged over the solid angle, cf. eq. (1). The largely negative values exhibited at radii 
larger than the average cage radius of 2.74 A indicate that a formally 4— charged Si cage 
would contain much more electron density at the outside as compared to the real VSi]*g 
system, while simultaneously there would be much less charge in the inside (positive regions 
in Fig. |3]). However, this does not simply indicate that a smaller formal charge transfer from 
metal to cage takes place. As illustrated in Fig. [3] also other superpositions of differently 
charged empty cages and cations do not represent the real electron density well. This holds in 
particular for the radial region between ~ 1 — 2 A, i.e. exactly the bonding region between 
central metal atom and cage. The metal-Si bonding is thus rather the result of a more 
complex hybridization than mere formal charge transfer. 

Equivalent results shown in Fig. Hlare also obtained for the "non-magic" TiSiJ^g and CrSi^^g 
clusters, which means that also there the real electron density of the endohedral cage cannot 
be fully rationalized in terms of a formal charge transfer. However, in all three dopant cases 
the true electron density outside the cage is best represented, i.e. the radial electron density 
difference distribution is closest to zero, for a charge combination of a positively charged Si 
cage and a neutral metal atom. This suggests that the different number of valence electrons 
in the three systems resides predominantly around the dopant. Figure O demonstrates that 
this is indeed the case. Depicted is the electron density difference between VSi^g and TiSi^g, 
as well as between VSi^g and CrSi^^g, which allows to locate the missing electron in TiSi^g 
and the excess electron of CrSij*g as compared to the "magic" VSi^g cluster, respectively. In 

9 



1 ■ 1 ' r i i ■ 1 1 1 i 1 -4-1 1 1 1 1 1 1 1 1 1 1 1 1 

0123456 012345 

Radial distance, A Radial distance, A 

FIG. 5. Electron density difference VSi^g - TiSiJj (left panel) and VSi^^g - CrSi^g (right panel). 
The radial electron density distribution, cf. eq. (1), as well as the 3D isosurface at 0.02 e/A^ in 
the inset demonstrate that the missing electron in the former and excess electron in the latter case 
are predominantly located around the central metal atom. 

both cases this is close to the central metal atom. 

A complementary view comes from the analysis of the projected DOS. For this Fig. |3p 
specifically compares the metal-projected DOS for the three doped cages. Interestingly, 
the metal contribution to the lower lying electronic shells up to the 2p shell is essentially 
the same in all cases. This is much different for the frontier shells 2d and l/i, which are 
mostly responsible for the bonding between cage and dopant. Here, there is a clear trend of 
increasing metal weight to the states when going from TiSi^g over VSi^'g to CrSij'g. If we add 
up these metal contributions over the occupied set of 2d and Ih states, we arrive at a total 
of 2.1 (Ti), 3.1 (V) and 3.7 (Cr) electrons in the three cases. Between Ti and V, as well 
as between V and Cr the metal dopant provides thus each time around one electron more 
to the hybridized states. The adapting degree of metal-Si hybridization hence compensates 
largely for the different total electron numbers. In other words, while from TiSi^g over 
VSi]*g to CrSi^^g there is each time one more valence electron in the topmost 2d and Ih 
shells, the number of electrons that is actually assigned to the Si cage remains essentially 
the same. The cage therefore effectively "sees" similar charge numbers, as the adaptive 
ability of the orbitals that are predominantly responsible for the metal-silicon bonding can 
accommodate for the charge excess or deficit. Intriguingly, it is, however, not just one state, 
e.g. intuitively the one with the changed occupation, that is responsible for this. Instead it 
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is the rehybridization of the entire set of 2d and Ih states, which effectively compensates for 
the " non-magicity" . This adaptive capabihty diminishes the role of electronic shell closure 
and is in our view the main reason that helps to stabilize the endohedral cage geometry also 
for TiSi^'g and CrSi^g despite their differing number of valence electrons. 



IV. CONCLUSIONS 



In summary our DFT-based unbiased configuration searches confirm the preceding inter- 
pretation of Lau et al."^^ that the Frank-Kasper polyhedron indeed represents the ground- 
state geometry for the series of doped TiSi^g, VSi^g and CrSi^g clusters. Endohedral doping 
can thus be used as avenue to stabilize cage-like Siig geometries. The electronic structure 
analysis demonstrates that all three systems are well described within the spherical potential 
model, i.e. the electronic manifold groups into states of defined radial and angular momen- 
tum number. Only the classic VSi^'g cluster achieves closure of the electronic 2d shell, while 
the varying number of valence electrons leads to an unoccupied 2d state in case of TiSi^^g 
and an occupied Ih state in case of CrSi^g. Shell closure is thus not a necessary condition 
for the stabilization of the cage-like geometry. 

We attribute this diminished role of shell closure for the stabilization to the adaptive 
capability of the metal-Si bonding, which is more the result of a complex hybridization than 
the mere formal charge transfer picture originally proposed in connection with the spherical 
potential model. This adaptive capability allows to locate the deficient electron in case 
of TiSi^^g, as well as the excess electron in case of CrSi^'g predominantly around the metal 
dopant. The effective charge assigned to the Si cage is then essentially the same in the three 
systems, i.e. the rehybridization of the 2d and Ih shells compensates for the "non-magicity". 
While electronic shell closure is still certainly a criterion for particularly enhanced stability, 
the flexibility of the metal-Si bond can help to stabilize also other cage-dopant combinations 
than predicted by this simple rule. This indicates the exciting prospect that a wider range of 
materials may eventually be cast into this useful geometry for cluster-assembled materials. 
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